We present a new solver for massively parallel simulations of fullythree-dimensional multiphase flows. The solver runs on a variety of computerarchitectures from laptops to supercomputers and on 65536 threads or more(limited only by the availability to us of more threads). The code is whollywritten by the authors in Fortran 2003 and uses a domain decomposition strategyfor parallelization with MPI. The fluid interface solver is based on a parallelimplementation of the LCRM hybrid Front Tracking/Level Set method designed tohandle highly deforming interfaces with complex topology changes. We discussthe implementation of this interface method and its particular suitability todistributed processing where all operations are carried out locally ondistributed subdomains. We have developed parallel GMRES and Multigriditerative solvers suited to the linear systems arising from the implicitsolution of the fluid velocities and pressure in the presence of strong densityand viscosity discontinuities across fluid phases. Particular attention isdrawn to the details and performance of the parallel Multigrid solver. The codeincludes modules for flow interaction with immersed solid objects, contact linedynamics, species and thermal transport with phase change. Here, however, wefocus on the simulation of the canonical problem of drop splash onto a liquidfilm and report on the parallel performance of the code on varying numbers ofthreads. The 3D simulations were run on mesh resolutions up to $1024^3$ withresults at the higher resolutions showing the fine details and features ofdroplet ejection, crown formation and rim instability observed under similarexperimental conditions. Keywords:
展开▼
